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1. Introduction 


Enzyme electrophoresis has been applied to the study of earthworm populations since the late 
1970s. Most studies have focused on allele and genotype variation in local populations. Surveys that 
have covered a wide geographic area are concerned with a few species. The introduced 
parthenogenetic, Octolasion tyrtaeum (SAVIGNY, 1826), has been surveyed in New York. 
Tennessee and North Carolina (JAENIKE er al., 1980; JAENIKE & SELANDER, 1979, 1985). There 
were 8 overall enzyme phenotypes. 2 of which predominated. The distribution of these types was 
attributed to random events, since there was no evident habitat preference. 

Dendrodrilus rubidus (SAVIGNY, 1826), a facultatively parthenogenetic species, also introdu- 
ced to North America was surveyed in the northeastern United States. The sample was made up of 
several enzyme phenotypes, but one-third of all individuals belonged to a single type. Some of the 
polyploid types may have originated in America (JAENIKE ef al., 1982). 

Eisenia nordenskioeldi (EISEN, 1879) is an endemic Palearctic species the distribution of which 
reaches the northern parts of the USSR. There are several parthenogenetic polyploid forms and a 
diploid amphimictic one (GRAFODATSKU ef al., 1982; PEREL’ & GRAFODATSKU, 1983 and 
BuLaTova et al., 1984). MALININA & PEREL’ (1984) studied enzyme variation in different 
polyploid forms of E. nordenskioeldi but did not assess their relative abundance. They found clonal 
co-occurrencies but did not record any clear ecological differences. 

COBOLLI SBORDONI et al. (1987) studied enzyme electrophoretic differences between triploid 
parthenogenetic and diploid amphimictic strains of the holarctic Allolobophora (= Aporrectodea) 
caliginosa (SAVIGNY, 1826) comjic\. Triploids in Italy and Algeria are genetically close to each 
other but deviate far from diploids in Italy and Great Britain. The parthenogenetic forms are 
probably autopolyploids. Studies have also been done on amphimictic species of the genera 
Aporrectodea (STILLE er al., 1980), Lumbricus (ROBOTTI, 1982a, b, 1984; ENCKELL et al., 1985a, 
b) and Eisenia (JAENIKE, 1982; ROBOTTI, 1982c and ØIEN & STENERSEN, 1984). In general these 
forms show a pattern of variation typical to bisexually reproducing animals. 

Dendrobaena octaedra (SAVIGNY, 1826) is an ecologically successful species with wide 
spectrum of habitats and found throughout Finland and northern Norway (TERHIVUO & VALO- 
VIRTA, 1978; TERHIVUO, 19882, 1989). The species prefers acid coniferous forest soils but is often 
found in wet habitats such as shores (BoucHÉ, 1972; TERHIVUO & VaLOVIRTA, 1978). It is 
considered to be polyploid (CASELLATO, 1987) and obligatorily parthenogenetic (OMODEO, 1955). 
There are hexaploid D. octaedra individuals in southern Finland (HONGELL & TERHIVUO, 1989). 

In the following we survey electrophoretic variation of D. octaedra in eastern Fennoscandia. 
We aim at assessing genetical affinities between populations and establishing whether the 
distribution of electrophoretically detectable enzyme phenotypes of the species is correlated to 
ecological and regional variables. 
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Fig. 1. Localities (1—85) for D. octaedra samples within the north-south (I—VI) and the east-west (VIII—IX) 
transects. The subareas (I—XI) are explained in the text. Vegetation zones are according to AHTI er al. (1968). 
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2. Materials and methods 


2.1. Sampling procedure, localities and biotopes 


The D. ocraedra individuals totalling 428 were collected from 85 localities along two transects, one running from 
north to south and another from cast to west. The numbers from 1 to 85 in fig. 1 show the order the localities sampled 
between July 2 and October 21, 1981. AII individuals from a locality were kept together in a glass bottle, provided 
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with an air supply and maintained in cold water containers. Upon arrival to the laboratory the worms were frozen and 
stored at about —24°C until assayed in October— November 1981. 

Samples were taken with 20—40 km intervals from each other along two transects. The north-south transect 
(localities 1—48) is about 1,200 km long. The east-west transect runs in southern Finland about 600 km from eastern 
parts of the country to the Áland Islands in the west (samples 49—85). The samples represented six site types called 
biotopes. The biotopes are as follows: 

1. Rich forests (RF) (28 localities/141 worms: loc. 2/7 worms, 4/5, 6/5, 13/5, 15/5, 16/5, 27/5, 31/5, 32/5, 36/2, 
40/5, 50/6. 53/5, 57/5. 58/5, 60/5, 61/6, 62/5, 65/5, 66/5, 67/5. 72/5, 74/5, 16/5, 78/5, 79/5, 82/5 and 83/5) refer to 
deciduous tree dominated forests or stands with abundant herbs, grasses and bushes. The soil is rich with mull and/or 
clay and receives annually a considerable addition of organic material such as falling leaves and decaying herbs and 
grasses. The pH and moisture content of the soil are higher than in coniferous forests. 

2. Myrtillus type coniferous heath forests (MT) (15 localities/80 worms: loc. 3/5 worms, 7/5, 49/5, 51/5, 54/5, 
55/5. 56/5, 59/5. 63/5, 64/5, 68/5. 69/5, 75/5, 77/5 and 85/10) have a podsolic soil profile with acid raw humus. The 
ground layer is covered with moss and Vaccinium myrtillus is the predominant shrub. 

3. Dry heath coniferous forests (DHF) (6 localities/30 worms: loc. 35/5 individuals, 38/5, 73/5, 80/5, 81/5 and 
84/5) of different types have well developed podsolic soil profile, low pH and low moisture content of the soil are 
characteristic to them. Different species of mosses and lichens dominate the field layer. Annual addition of organic 
matter to the soil is low; much of it consists of pine needles. 

The sequence RF-MT-DHF represents a descending gradient of luxuriance of vegetation, soil moisture, pH and 
quality of organic material that earthworms may utilize. The Lumbricidae species assemblages in these forests differ 
widely from each other and due to higher density of individuals and biomass and as argued in TERHIVUO (198823) 
between species competition may be higher in RF than either in MT and DHF. 


Table 1. Grouping of localities (n = 85) and D. octaedra individuals (n = 428) sampled along the north-south (N-S) 
and east-west (E-W) transects. 


Sub- Trans- Serial No of No of 
area ect loc. no loc. worms 
I N-S 40—46 7 31 
i N-S 37—39, 47 4 19 
IH N-S 30—36. 48 8 37 
IV N-S 22—29 8 42 
V N-S 12-21 10 50 
VI N-S 2—11 10 52 
Vu E-W 54-61 8 41 
VIII E-W 50—53, 62 5 26 
IX E-W 1, 49, 63—68, ll 60 
75—76, 85 
X E-W 69—74 6 30 
XI E-W 77—84 8 40 


Note: The localities are indicated in fig, 1. For further explanation see the text. 


4. Shores (Sh) (15 localities/70 worms: loc. 1/5 worms, 5/5, 18/5, 34/5. 37/5, 39/5, 41/5, 43/2, 44/4, 45/5, 46/5, 
47/4, 48/5. 52/5 and 70/5) refer to both alluvial sea shores by the Gulf of Finland and to river banks and lake shores in 
other parts of the study area. The sea water is brackish with 0.4—0.5% salt. The soil is wet. occasionally waterlogged 
but often rich in plant debris brought by waves or flooding. 

5. Meadows (Me) (19 localities/96 worms: loc. 8/5 worms, 9/5, 10/5, 11/5, 12/5, 14/5, 17/5, 19/5, 20/4, 22/6, 
23/5, 24/5, 25/5, 26/5, 28/6, 29/5, 30/4, 33/6 and 42/S) are characterized by the abundance of grasses and herbs. Trees 
and bushes are in general absent. The organic soil layer is thick, rather moist and rich with mull and/or clay and there is 
often silt underneath. In Finland meadows host higher numbers of earthworm individuals per square metre than any 
other biotope (TERHIVUO, 1989). 

6. Waste or fallow soils (WS) (2 localities/1 1 worms: loc. 21/6 worms and 71/5) are sites with high human impact 
upon both the content and structure of the soil. These sites harbour ruderal plants. The soil is either fertile with mull 
(locality 71) or poor with a thin layer of plant debris and silt underneath (loc. 21). 
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Fig. 2. Relative frequencies (% ) of sample localities representing 6 types of biotopes in the subareas along the N-S 
and E-W transects. Symbols for the biotopes: RF— rich deciduous forests, MT = Myrtillus type coniferous moist 
heath forests, DHF = dry heath coniferous forests, Sh — shores, ME= meadows and WS = waste soils near human 
settlements. 


For further information of the biotopes see Aario (1960) and TERHIVUO (1989). 

To compare parts of the transects the sample localities were poststratified into stratum called subareas (fig. 1, 
table 1). The subarea stratum consists of samples in 100 km X 100 km squares of the Finnish uniform grid system 
(HEIKINHEIMO & RAATIKAINEN, 1971). 

Fig. 2 shows the proportions of the different biotopes within the subareas. Clearly the biotopes are not equally 
represented along the transects. To assess within and between biotope variabilities in different parts of the two 
transects the following categories for the samples are used: the northern part of the N-S transect is called Lapland 
(subareas 1— III), central Finland comprises subareas IV—V. together with Lapland they constitute a transect unit 
called North. The E-W transect is made up of three parts. East (subareas VII—VIII), West (subareas IX—X) and Aland 
(subarea XI). Together with subarea VI they form another transect unit called South. The Lappish localities are in 
general within the North Boreal vegetation zone and the ones in central Finland are located to the Middle Boreal zone 
(fig. 1). The southern localities are divided between the South Boreal and Hemiboreal zones. 


2.2. Electrophoresis and enzyme assay methods 


A set of 4—7 anteriormost segments were cut off from each worm and homogenized. The homogenate was 
assayed for the following enzymes: esterases (Esr-/ and Esr-2), glutamate-oxaloacetate transaminase (Got), 
glyceraldehyde-3-phosphate dehydrogenase (G-3-pdh), isocitrate dehydrogenase (/dh), leucine aminopeptidase 
(Lap), malate dehydrogenase (Mdh), phosphoglucomutase (Pgm) and phosphoglucose isomerase (Pgi) using 
electrophoresis and assaying methods described by SAURA er al. (1979). Esr-7 and G-3-pdh were invariable in the 
entire material. The bands of Got proved to be poorly reproducible so that it is not used in defining the overall enzyme 
phenotype, either. 

The enzyme systems adopted are the same used in studies of other earthworm species (e.g. JAENIKE er al., 1980; 
STILLE et al., 1980). It appears that the numbers of allozymes within a given enzyme system are different between 
species. 

The most frequent band of an enzyme in the assayed gels is given the index 100 and the other bands are identified 
and correspondingly indexed according to their millimetric difference in anodal migration from that band. The band 
patterns observed are indicated in table 2. 
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Table 2. Enzyme systems surveyed and band patterns of the enzyme variants recorded in 428 D. octaedra worms 
from 85 localities in Finland and N Norway. The numbers of populations (= localities) for the variants are also 
indicated. 


Enzymes and serial Band pattern No of Of all 
number of variants loc. loc. 
TUE C ee a n [%] 

Mdh l 1.00/1.05 80 94.1 

2 0.97/1.00 59 69.4 

Est-2 3 1.00 20 23.5 

4 1.00/1.05 57 67.1 

5 0.95/1.00/1.05 54 63.5 

6 1.00/1.03/1.05 16 18.8 

7 0.95 5 5.9 

8 0.95/0.96/0.98 17 20.0 

9 0.95/1.02/1.07 7 8.2 

10 0.94/1.05 Il 12.9 

1 0.85/0.90/0.95/1.00 14 16.5 

12 0.87/0.97/1.05 11 12.9 

13 0.87/1.07 l 1:2. 

14 1.05 3 3.5 

15 0.94/0.98/1.00 8 9.4 

Pgm 16 1.00/1.05 85 100.0 

17 0.95/1.05/1.15 5 5.9 

18 1.00/1.05/1.15 15 17.6 

idh 19 1.00 84 98.9 

20 1.00/1.07 31 36.5 

Lap 21 1.00 66 77.6 

22 1.00/1.15 43 50.6 

23 0.85/1.00/1.10/1.15 4 4.7 

24 0.98/1.00/1.15 l 1.2 

25 0.96/1.00 6 7.1 

26 1.00/1.03 1 1.2 

27 1.00/1.02 3 3:5 

Pgi 28 1.00 79 92.2 

29 1.00/1.03/1.06 E 4.7 

30 1.00/1.04/1.08 6 7.1 

31 0.94/0.97/1.00 30 35.3 

32 0.89/0.95/1.00 10 11.8 


Note: The numbers of populations (= localities) for the variants are also indicated. The band patterns are explained in 
the text. 


2.3. Parameters calculated 


Clone pool diversities of D. octaedra within and between the biotopes and along the transects were studied by the 
rarefaction method. This method avoids the difficulty of dissimilar sample size by predicting the number of clones for 
samples smaller than the ones in the original data set. The formula is as follows: 


S (Nm) 
E(S,) = * |1- 
i=! (=) 
\n 


in which E(S,,) is the expected number of clones in a sample of n individuals derived from the total data set, S is the total 
number of clones, N denotes the number of all the individuals, N; is the number of individuals in the ith clone and n is a 
random sample of n individuals (SANDERS, 1968; SIMBERLOFF, 1978, 1979). 
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JANSON & VEGELIUS (1981) proposed a set of six criteria that an efficient similarity measure (r) must fulfil. They 
found that out of 20 similarity measures that they studied, three conformed to the prerequisites that they proposed. 
They are: 

(1) The formula of OcH1at (1957) 


A 
VA+BxVAFC 


where A = the number of clones shared by areas v, and v2, B = the number of clones present in area v; but not in v2, 
C = the number of clones present in area vz but not in v}. 
(2) The formula of Dice (1945) 


ip = 


2A 
Tp = —————— 
2A+B+C 

where A, B and C are as above. 
(3) The formula of Jaccard (1908) 

A 
Tj = = 
A+B+C 


where A, B and C are as above. 

These three formulas were applied to assess clone poo! similarities between different transects units and different 
biotopes. 

Symbols for significances of the statistical tests applied are indicated in the legend to table 4. 


Table 3. Symbols for overall enzyme phenotypes of the D. octaedra clones with four or more individuals recorded in 
a sample of 428 individuals from different parts of Finland and N Norway. 


Enzyme Symbols for clones and combinations of enzyme variants 


yen (à (B (CO (D (E (b) (O (HG (D QO (D (D (Mw 
Mdh 1 1 1 1 l 2 23 23 2 1 1 1 1 
Est-2 5 4 4 3 5 5 4 3 4 3 6 4 ll 
Pgm 16 16 16 16 16 16 16 16 16 16 16 16 16 
Idh 19 19 19 19 19 19 19 19 19 19 19 20 19 
Lap 21 21 22 22 21 21 22 21 21 21 21 22 21 
Pgi 28 28 28 28 28 28 28 28 28 28 28 28 28 
No of 

ind. 40 34 25 15 15 13 10 10 8 7 7 7 6 
No fo 

loc. 21 20 12 9 11 11 5 6 7 7 7 4 5 
Enzyme 

system NW © ® @ ® G (D wo wv o WM o 
Mdh 1 1 ] 1 2 1 1 1 2 1 2 2 
Est-2 3 1 0 8 I 4 8 2 S 4 3 4 
Pgm 6 6 6 6 6 6 6 6 6 (18 16 16 
Idh 9 9 19 #19 19 20 19 19 19 20 19 20 
Lap 22 21 21 21 21 21 22 22 22 21 22 21 
Pei 2 31 28 28 28 028 028 2% 2 2 2 28 
No of 

ind. 6 6 5 5 5 4 4 4 4 4 4 4 
No of 

loc. 1 3 35 4 3 232 4 23 ^4 2 1 4 
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3. Results 
3.1. Origin, range and ecology of D. octaedra clones 


3.1.1. Origin and extent of variation in enzyme systems 


Six loci (Est-2, Idh, Lap, Mdh, Pgm and Pgi) proved to be variable. The extent of variation is 
shown in table 2. A total of 147 different overall enzyme phenotypes (or clones for short) were 
recognized in the material. Out of the 147 clones, 80 were unique, i.e. found only once. The 18 
most common clones make up 52.3% of all the worms. These clones are identified with letters from 
A to R and their allele constitutions are shown in table 3 together with those of clones S—Z each with 
4 worms. 

The D. octaedra clones may be derived from each other through a single mutation step (i.e. 
through an addition or deletion of an allele in any enzyme system) or through some recombination. 
Because the worms are also polyploid attempts to establish the nature of the genetic change leading 
to a different enzyme phenotype can not be considered worthwhile. Anyway, an addition of a novel 
allele should be visible, as well as the loss of an allele in a situation, where there are no other alleles 
of identical allozyme product. Polyploid parthenogenetic populations may also harbor null alleles, 
which make the dose of an allele (e.g. five copies of allele A, one copy of allele B, as judged from 
the enzyme phenotype) difficult to assess. We have, therefore. given up efforts to indicate any 
dosage effect which might underly the enzyme variants. 

D. octaedra reproduces through apomictic parthenogenesis while other parthenogenetic 
earthworms are automictic (MULDAL, 1952; OMopEo 1952a, 1955). The genotype is thus 
transmitted intact to the offspring without any recombination of alleles. One can, therefore, 
consider different loci only in the context of the genotype as a whole. To assess the contribution of 


Table 4. The statistical (4?) significances of polymorphic D. octaedra enzyme marker systems revealed between and 
within the transects and biotopes. 


Polymorphic enzyme marker group 


Mdh Est-2 Pgm Idh Lap Pgi 

Transects 

1. Between * ee = kkt dide o 
transects 

2. Within N-S doe = * * Es E 

3. Within E-W = = = o ek * 

Biotopes 

4. Forests/open biotopes : o = = +x ** m 
(totals of variants) i 

5. Variants between - = - xc ** ** 
forest types 

6. Totals of variants - * = = ek ee 
between forests of 
north and south 

7. Variants between types sxe o - [o] = * 
of open biotopes 

8. Totals of variants sk = = * = = 


between open biotopes 
of north and south 


Note: Symbols: *** = P=0.001, ** = 0.01=P>0.001, * = 0.05=P>0.01 and O = 0.10=P>0.05. 

In the within transect comparisons the samples are pooled according to subareas and x?-test was applied to the 
corresponding frequencies of the most common enzyme variants over the transects. In the between transect 
comparisons all the samples are correspondingly pooled. As for biotope comparisons the samples are pooled within the 
transect units (north-south, see chapter 2.1. and fig. 1) or between the two transects. 
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Fig. 3. Distribution patterns for D. octaedra clones with = 5 individuals: (1) This pattern involves clones A, B, C, F, 
1, K&L. (2) Clones D, H and J. (3) Clones E, M, P&Q. (4) Clones G, N, O&R. For further information see the text. 
For enzyme variant combinations characterizing the clones see table 3. 


individual enzyme systems to the overall enzyme phenotypes the enzyme variant markers 
responsible for the clonal diversity in different parts of the transects and between the biotopes were 
evaluated. All the polymorphic enzymes are dissimilar in this respect (table 4). 
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3.1.2. Distributions and ecological divergence of common clones 


The geographical distributions of the common clones represented by at least five individuals 
show dissimilar patterns: (1) some extend along both N-S and E-W transects, (2) some cover the 
E-W transect as well as the southern and central parts of the N-S transect, (3) some extend along the 
E-W transect to southern parts of the N-S transect and (4) some can be found only in the central part 
of the N-S transect (fig. 3). The most abundant clone ist type A, which was represented by 
40 individuals from 21 localities all over the study area. Many of the distributions for the common 
clones may be explained on the basis of random events, but, considering e.g. clone E recorded in 
11 localities all lying within Hemiboreal and South Boreal vegetation zones, this may not be the 
case. There are 45 localities out of 85 within these two zones and the probability to have within them 
11 localities of clone E by change is thus (45/85)! , i.e. P« 0.001. 

Of all the lumbricids inhabiting Fennoscandia D. octaedra occupies the widest spectrum of 
habitats (TERHIVUO & VALOVIRTA, 1978). The adaptability of a parthenogenetic taxon is the total of 
the adaptabilities of its component clones. We approached the biotope preferences by considering 
the presence of the most abundant clones in each biotope in different parts of the transects. 

Table 5 shows how the most abundant clones representing different distribution patterns in fig. 3 
occur in the six biotopes and in the two transects units. Most of these seem to be broad-niched but 
they are unequally distributed between the biotopes. For instance, the most common clones A, B 


Table 5. Habitat distributions of some common D. octaedra clones in different transect units. Symbols for clones 
correspond to those in table 3. 


Biotope Clones Most 


and ee frequent 
transect unit a B e D E G clone 
Forests (all) 21 28 10 11 12 3 B 
north 3 9 8 3 -— 3 B 
south 18 19 2 8 12 = B 
Rich forests 11 16 9 7 8 2 B 
north 3 4 8 3 - 2 c 
south 8 12 l E 8 - B 
Conif. forests 10 6 — 4 4 1 A 
north - 5 = = = I B 
south 10 7 1 4 4 - A 
Open biotopes (all) 9 6 15 4 3 7 A 
north 10 3 1 2 = - B 
south 9 3 4 2 3 7 A 
Shores 13 3 8 2 l - A 
north 9 2 7 = = - A 
south 4 1 l 2 l - A 
Meadows 5 3 7 -— - 1 C/G 
north - 1 4 = = F G 
south 5 2 3 = = = A 
Waste soils 1 = = — 2 = E 
north = = = = — 2€ = 
south l = = = 2 = E 
Total 40 34 25 15 15 10 A 
Forests (%) 53 82 40 73 80 30 B 
Open biot. (%) 47 18 60 27 20 70 A 
Distrib. type l 1 i 2 3 4 — 


Note: For the areal distribution types of the clones A—G see the text and fig. 3. 
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Fig. 4. Changes in relative frequencies (96) of the enzyme variants present in the most frequent D. octaedra clone A 
and in groups of clones that deviate in one to 4 enzyme variants from the corresponding combination present in thc 
clone A. In clone A the enzyme variant compilation ist Mdh (= variant 1)-Est-2 (5)-Pgm(16)-Idh(19)-Lap(21)- 
Pgi(28) thus equalling 100% in the diagram. For further explanation see the text. 
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Fig. 5. Mean enzyme variant differences (X + — S.E.) in D. octaedra individuals from the subareas as measured by 
the numbers of the enzyme variants that they differ from those present in clone A. In the diagram the individuals of 


clone A are included and considered equal to 0. If these are excluded highly significant statistical correlation 
(r=0.933***, df=9) to the diagram shown here still exist. See also the text. 
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Fig. 6. Relative frequencies (% ) for the individuals of the most abundant D. octaedra clone A as well as for those 


deviating in one (G1) or 2 (G2) enzyme variants from it withing the subareas. Broken lines indicate the mean for each 
clone group as counted over the means of the subareas. 


and C which were recorded all over the study area are unequally distributed between forests and open 
biotopes (y^ = 12.136**, df = 2). 


3.1.3. Distributions and biotope divergence of clone groups 


Clone A is the most frequent overall phenotype. It was found in all biotopes and different parts of 
the transects and may be taken to represent a successful general-purpose type. If clone A is taken as a 
kind of East Fennoscandian D. octaedra norm, from which the other types can be derived, we may 
consider the distribution and biotope preferences of the other clones in relation to A. The grouping 
allows, however, another way surveying clones in relation to regional and biotope variables. In this 
study the minimum number of enzyme variation steps needed to derive a clone from A is between one 
and four. This, of course, does not mean that the other clones necessarily are derived from it. 

The first group comprises 18 clones and 121 worms from 57 localities; the second 52 clones/142 
worms/69 localities, the third 60 clones/104 worms/49 localities and the fourth 16 clones/21 worms 
from 16 localities. The mean number of worms per clone in these groups is 6.7, 2.7, 1.7 and 1.3, 
respectively. In other words, in the first group about every seventh worm represents a novel genotype 
butin the fourth group about all worms are dissimilar in this respect. Most clones in group 1 differ from 
A in the Est-2 locus, while in the other groups the difference is due to other loci as well (fig. 4). 

Fig. 5 shows the proportionate genetic difference of individuals from clone A in each subarea. 
Worms of clone A are included but they are given the value zero. In all subareas of the N-S transect the 
means are higher than along the E-W transect. The highest difference is found in subarea IV. The 
result can be explained either through differences in adaptation or migration history or both. 
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Fig. 7. Relative frequencies (% ) for the individuals of the D. octaedra clone groups that deviate in 3 (G3) or 4 (G4) 


enzyme variants from those of the most abundant clone A in each subarea. Broken lines indicate the mean for each 
clone group as counted over the means of the subareas. 


Fig. 6 shows the relative frequencies of worms of clone A and clone groups 1 and 2 in different 
subareas. Clone A is distributed throughout the subareas regardless of their associated biotope. 
Their relative frequency never exceeds 20%. Group 1 dominates along the E-W transect and 
group 2 along the N-S transect. While the relative proportions of clones fluctuate the difference 
between two clone groups is uniform. 

Fig. 7 shows the relative frequencies of worms of groups 3 and 4 along the transects. They are 
less common than worms of the preceding types, in particular along the E-W transect. In fact, one 
may note the relative rarity of group 3 and the virtual absence of group 4 in the south. 

A comparison of the frequencies of the four clone groups and that of clone A in the two transects 
shows a highly significant (y? = 29.261***, df = 4) difference. In other words, they are unevenly 
distributed. According to figs. 6 and 7 worms of the Hemiboreal or South Boreal vegetation zone 
have diverged less from clone A than the northern worms. 

In forests D. octaedra worms of clone groups 1—4 differ from clone A in about 1.89 enzyme 
steps but in open biotopes by 2.32 enzyme steps (table 6). Worms of different clone groups are not 
homogeneously spread over the two main biotope categories (y^ = 24.409***, df = 3). In addition, 
in forests the difference between the north and the south is not significant for clone groups 1, 2 and 
the combined clone group 3—4 (3? = 2.307, ns, df = 2) but there is a difference in open biotopes 
(y) = 15.360***, df = 2, groups 3 and 4 combined). 

The distributions of clone groups between the three forest biotopes are equal over the transect 
units (X? — 0.529, ns, df — 4, groups 3 and 4 pooled). Considering meadows and shores slight 
heterogeneity is indicated (y? = 4.9343°, groups 3 and 4 pooled). Worms of the clone group 4 seem 
to account for this. It may be noted that they are absent from the south, except from the Áland 
islands. 


124  Pedobiologia 34 (1990) 2 


Table 6. Areal and biotope distributions of D. octaedra clone groups deviating in one to 4 enzyme types from those in 
clone A. 


Biotope and No. of | Clone group Difference in 
transect unit indiv. > «= — 9 24a -  — enzyme steps 
from clone A 


n [%] m [99] n [%] m [5] RSE. 

Forests je. 230 9| 396 80 348 53 230 6 2.6 1.89 + 0.06 

(grand 1. RF 130 51 392 47 362 27 208 5 3.8 1.89 + 0.07 

total) 2. MT 70 27 386 24 343 19 241 - —  18950.0 

3.DHF 30 13 433 9 300 7 233 1 3.3 1.87 € 0.16 

north I. RF 34 l0 294 13 382 9 265 2 59 2.09 +0.15 
LMT = = = SS = = = + uS = 

3. DHF 10 6 600 2 200 2 200 - — 160-027 

total 44 16 364 15 341 I 250 2 46 1.89+0.14 

south l. RF 96 41 42: 34 354 18 188 3 3.1 1.82 0.09 

2. MT 70 27 2386 24 343 19 231 - — 1890.10 

3.DHF 20 7 350 7 350 5 250 | 50 2.00 + 0.21 

total 186 75 403 65 35.0 42 226 4 2.20 1.87 + 0.06 

Open 4-6. 158 30. 19.0 62 39.2 51 323 15 9.5 2.32 + 0.07 

biotopes 4. Sh 5 12 2LI 29 309 15 263 | L8 2.09 + 0.10 

(grand 5. Me 91 16 20.9 33 363 30 330 12 132 2.42+0.10 

total) 6. WS 10 2 39 = «= 6 60.0 2 20.0 2.80+0.33 

north 4. Sh 4l 5 122 2 S512 4 341 | 24 227 £0. 

5. Me 75 ll 147 26 347 28 373 10 133 2490.1 

6. WS ^ = - - = 4 667 2 333 3335021 

total 122 16 13.1 47 383 46 37.7 13 107 2.46 * 0.08 

south 4. Sh 16 7 438 8 500 | 63 - = 16340.16 

5. Me 16 5 313 7 438 2 125 2 125 2.06 +0.25 

6. WS 4 2 500 = - 2 500 — — 2.00+0.58 

total 36 ij4 389 15 47 5 idi» 5 56 1.86 +0.14 

Total (1—6.) 383 121 31.2 142 36.6 104 26.8 21 5.4 2.06 + 0.05 


Note: For further explanation see the text. 


3.1.4. Unique clones 


The clones represented by a single worm are unique. They total 80 and make up 18.7% of all the 
worms but 54.4% of all the clones. Unique clones may represent marginal populations of more 
widespread clones, novel mutants that have recently arisen in previously existing clones or products 
of recombination events. Mutants may occur in populations with about equal frequency and they 
should be evenly distributed; recombination events might be expected to occur only in some 
populations prone to them and thus show an areal impact. 

The frequencies of unique clones were not significantly different between subareas (y? = 
12.062, ns, df = 9), subareas 8 and 9 pooled, (table 7) but were between the biotopes (i? = 
12.179*, df = 5) (table 8). There is also a slightly significant difference between forests and open 
biotopes (4^ = 2.977°, df = 1) as well as between open biotopes in the south and the north (^ = 
3.158, df = 1). In the forests unique clones are, however, evenly distributed. 

The number of unique clones is by far highest in the clone group 4, which deviates genetically 
most from clone A. Unique clones make up in group 4 as much as 57.1 % of individuals and 75 % of 
clones. Unique clones are few in clone group 1, where they make up only 4.1% of individuals and 
27.8% clones. while groups 2 and 3 are intermediate between the two extremes. Those unique 
clones that deviate most from clone A are most frequently recorded in central parts of Finland. 
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Table 7. Proportions of D. octaedra individuals representing rare (= unique) clones in different subareas. 


Subarea Total No of 

of unique 

indiv. clones 

* n [96] 
I 31 6 19.4 
Il 19 5 26.3 
lll 37 7 18.9 
IV 42 6 14.3 
V 50 15 30.0 
VI 52 5 9.6 
VII 41 10 24.4 
Vill 26 4 15.4 
IX 60 8 13.3 
X 30 8 26.7 
IX 40 6 15.0 
Total 428 80 18.7 


Table 8. Frequencies of unique D. octaedra clones in different biotopes and transect units. 


Biotope No. of No. of Total 
unique worms in number 
clones other clones of worms 
n [95] n [99] s 

Forests (total) 40 15.9 211 84.1 251 

north 5 11.9 42 88.1 47 
south 35 25.0 169 75.0 204 
Rich forests 22 15.6 119 84.4 141 
north 4 10.8 33 89.2 37 
south 18 17.3 86 82.7 104 
MT forests (= south) 12 15.0 68 85.0 80 
Dry heath forests 6 20.0 24 80.0 30 
north 1 10.0 9 90.0 10 
south 5 25.0 15 75.0 20 
Open habitats (total) 40 22.6 137 77.4 177 
north 34 25.8 98 74.2 132 
south 6 13.3 39 86.7 45 
Shores 12 17.1 58 82.9 70 
north 9 18.0 41 82.0 50 
south 3 15.0 17 85.0 20 
Meadows 22 22.9 74 77.1 96 
north 21 21.9 55 78.1 76 
south 1 5.0 19 95.0 20 
Waste soils 6 54.5 5 45.5 11 
north 4 66.7 2 33.3 6 
south 2 40.0 3 60.0 5 
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Fig. 8. Frequency distributions of electrophoretically determined D. octaedra clones in samples ( — localities) with 
unequal numbers of individuals. 


The proportion of unique clones was more variable between open biotopes (range 5— 66%) than 
between the forest types, where the frequency of unique clones was rather even (range 10—25%) 
(table 8). 


3.2. Diversity of D. octaedra clones 


3.2.1. Within sample diversity of clones 


Fig. 8 shows the within sample variability of clones in relation to sample size. Within three 
samples all five worms were identical with each other. In all other samples at least one worm 
differed from the others. Occurrence of worms with similar allozymic combinations evidences 
parthenogenetic reproduction within populations. 
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3.2.2. Diversity of clones along the transects 


Clonal diversity was lowest in northern parts of the study area (subareas I—IV) but no decrease 
was observed along the E-W transect (table 9). In central Finland (subareas V — VI) the diversity of 
clones was, however, rather high. This may be due to post-glacial migration of D. octaedra 
populations from N Sweden. One can expect to find about 25—29 clones in a sample of 40 worms 


Table 9. Predicted numbers E(S) of D. octaedra clones calculated for rarefied samples of 30 individuals and 
40 individuals from different parts of the transects. 


Subarea 30 ind./sample 40 ind./sample Transect 
pup E(S) Range E(S) Range 

I-II 19.8 18.2—21.3 23.7 22.4—25.0 N-S 

III -IV 20.4 18.6—22.3 24.8 22.8—26.7 N-S 
V-VI 22.4 20.4—24.4 27.9 25.7—30.2 N-S 
VII- VIII 20.8 18.9—22.7 25.9 24.0-27.8 E-W 
IX-X 23.0 21.1—24.9 28.9 26.8—31.1 E-W 
XI 20.6 19.3—-21.9 25.0 (total catch) E-W 


Note: Range = [E(S) — S.D.] — [E(S) + S.D.]. Standard deviation = S.D. For further explanation see text. 


Table 10. Predicted numbers E(S) of D. octaedra clones calculated for rarefied samples of 10, 20 and 30 individuals 
from different biotopes and transect units. 


Biotope and Predicted number fo clones 
— 10 ind./sample 20 ind./sample 30 ind./sample 
E(S) € S.D. E(S) + S.D. E(S) + S.D. 

Forests: 

Rich forests (total) 8.8 + 1.0 15:9 :& 1.7 31,9 x 2.1 
north TEST 11: x 4:2 13.6 + 1.0 
south 8.9 + 1.0 16.0 + 1.6 22.2:F 2.0: 
-east 8.8 + 0.9 15:8 £13 214: 1.2 
-west 8.7 € 1.0 15.8 € 1.2 22 
-Aland 8.7 + 0.9 15 - 

MT-forests 8.8 + 1.0 15.9 + 1.6 21.9 + 1.9 

(total — south) 

-east 8.4 + 1.0 14.5 + 1.0 = 
-west 8.7 € 0.9 15.5 = 1:3 214 3E 12] 
-Aland 9 za 

Dry heath 7.8.:E 1.1 13:0 = 1:2 17 

forests (total) 
north 5 = = 
south 7.9+1.0 13 = 

Open biotopes: 

Shores (total) 8.2 d 1.1 14:335 1:7 19.3 1.9 
north 8.0 + 1.1 13.5 + 1.6 18.0 + 1.6 
south 7.9+ 1.0 13 = 

Meadows (total) 9.2+0.8 16.9 + 1.4 23.6 + 1.9 
north 9.2+ 0.8 17.1 1.4 24.0 + 1.8 
south 7.6 + 1.0 12 = 

Waste soils (total) 9 - = 


Note: The figures without standard deviation (S.D.) indicate the actual numbers of clones found. 
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taken in southern and central parts (subarea V — XI) of Finland (table 9) but in northern Finland and 
Norway (subareas I—IV) the expectation is 23—25 clones. 


3.2.3. Diversity of clones in different biotopes 


According to fig. 2 samples from different biotopes are unequally distributed along the two 
transects. One is thus confronted with two unequal data sets. To overcome this difficulty the 
rarefaction method is applied to samples pooled according to biotope and units of the transects. 
Table 10 shows the expected total number of clones for each biotope based on the information 
gained from all the samples of the biotope. i 

The data for the areal units within the biotope are based on corresponding sample sets, but as this 
information has already once been used in measuring the totals table 10 merely warrants 
comparisons between the totals or between parts of transects, respectively. Results from the 
rarefaction suggests that in a sample of 30 individuals about 22 different clones can be found in RF 
and MT forests but only about 17 in DHF (table 10). Within the RF and DHF forests the diversity 
decreases towards the north, at least to some extent. In general, DHF forests host fewer clones than 
RF and MT forests do. 

In meadow habitat the number of clones is lower in the south in comparison with the north. 
Meadows seem to have higher number of clones than shores in any part of the transects. 

In conclusion, the number of clones decreases towards the north in RF and DHF forests but 
increases in meadows. Different habitats may host a different array of clones. The clone pool of 
D. octaedra thus varies according to geographical and environmental factors. 


3.3. Similarities in clone associations 
3.3.1. Areal similarities 


The diversity measures used in the two preceding sections are based on the total number of 
clones. Diversities do not, however, indicate the proportion of clones that are shared between areas 
or biotopes in different parts of the two transects. Similarity is just as interesting as difference; it 
may provide evidence to aiding the interpretation of the adaptation and migration history of the 
species. 

Similarity measures applied in ecology most often involve cases of coexistence of two species 
over different localities or, alternatively, similarities measured between two localities over different 
species. In this study, the latter approach was applied to clones instead of species. 

The formulas of OcHIAI, DICE and JACCARD (see section 2.3.) were applied to the data on clonal 
occurrences within each two subareas in pairwise comparison. The three indices are highly 
significantly (P<0.001) correlated to each other, as the following tabulation shows (df = 43 in 
each comparison): 


OcHIAI JACCARD 
(To) (ry) 
DICE (rp) 1.000 0.998 
JACCARD (rj) 0.995 — 


As this is the case, only one formula needs to be utilized in presenting the data. Fig. 9 presents 
the results according to the index of Ocuiar. In fig. 9 the index measures correspond to the 
polygonal matrix of similarity values between transect units compared in pairs. The similarities are, 
on average, very low and do not exceed 0.5 in any comparison. Units of the E-W transect show 
higher affinities to each other than to central or northern parts of the N-S transect. This supports the 
contention that the clone pools along the transects are different. 
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Fig. 9. Similarities in D. octaedra clone pools between the subareas as shown by the index of OcH1A1 (1957). See also 
the text. 
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Fig. 10. Similarities in D. octaedra clone pools between the biotopes sampled in northern and southern parts of the 
transects as measured by the index of OcHi1A1 (1957). See also the text. 
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3.3.2. Similarities in clone associations between biotopes 


Similarities between the clone pools from different biotopes are low (fig. 10). In the south rich 
forests and MT forests as well as meadows and shores share the highest number of clones. In the 
north, however, rich forests show affinities to most other habitats and meadows share clones with 
shores. 

The low proportion of clones in common between biotopes in the north and the south implies 
that a majority of clones are adapted to inhabit a given biotope and/or a geographic area. 

The similarities observed here are mainly due to clones with a wide geographic range and a 
broad spectrum of habitats, that is, clones that can be designated as eurytopic. The majority of 
individuals seem to belong to stenotopic clones. Taken together, this is an indication of the 
evolutionary potential of a parthenogenetic species living in the northern margin of its geographic 
distribution. ` 


4. Discussion 
4.1. Genetic divergence of clones 


In most cases the common clones can be derived from clone A through single mutation steps. 
Fig. 11 shows an attempt to tie the relationships of these clones to each other and to clone A as well 
as the number of worms found. It has been customary to give such relationships in the form of a 
dendrogram but we here attempt to show what kind of changes can be involved. Evidently all other 
loci but Pgm have mutated to produce the less common types shown here. As in the case of 
dendrograms, this figure is open to other interpretations. For example, clone N may originate from 
clone A by way of D or J. Clone C can be derived from A through B or D etc. 

Omopeo (1953) argued that the apomictic polyploid D. octaedra reproduces exclusively 
through parthenogenesis. Chromosomal evidence by HONGELL & TERHIVUO (1989) implies high 
degree of allozyme polyploidy in eastern Fennoscandian populations and supports this contention. 
The pattern of variation observed in eastern Fennoscandian populations indicate that D. octaedra 
clones represent a monophyletic lineage. The species is morphologically highly variable (OMODEO, 
1952b; Gates, 1974) but attempts to correlate genotype differences with morphological variability 
were not successful (TERHIVUO, 1988a). 


4.2. Enzyme markers in relation to clone pools 


CHRISTENSEN (1980), JAENIKE (1982), ØIEN & STENERSEN (1984) and MALININA & PEREL’ 
(1984) showed that as few as one or two enzyme markers are useful characteristics in separating 
groups of parthenogenetic earthworm clones or species from each other. 

Fig. 12 shows a hypothetical case of a polymorphic clonal population in relation to a gradient. 
We assume that the clones respond dissimilarly to the gradient so that some become more common 
while the others grow scarce. Each clone is characterized by marker systems, some of which 
indicate no response to the gradient (A1— A3), others show cleat-cut differences in different parts of 
gradient (B1—B2; C1 and C2) and the rest exhibit trends rather than discontinuities along the 
gradient. Fig. 12A shows a case where no response to the gradient can be seen. In fact the clones 
respond but the markers fail to show it. Fig. 12B illustrates a situation of graded response to the 
gradient, similar to that reported in Lumbricillus lineatus (CHRISTENSEN et al., 1976, 1977; 
CHRISTENSEN, 1980). In fig. 12C different clones exploit resources along the gradient unequally. 
One now observes trends rather than discrete differences between the clones. This case comes 
closest to our interpretation of the data on D. octaedra. 

JAENIKE & SELANDER (1985) surveyed clones of the parthenogenetic earthworm Octolasion 
tyrtaeum without finding evidence for ecological differentiation between clones. They attributed 
the clonal variability to stochastic factors such as clonal drift. This does not, however, explain the 
fact that they found a regular coexistence of predominant clones in an extensive material. 
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Biotype A 
V (40 indiv.) 


Fig. 11. Numbers of individuals and the enzyme systems involved in the genotype dissimilarities between some 
abundant D. octaedra clones as compared to those for the most abundant clone A. The symbols denote the 
polymorphic enzymes but not any individual variant responsible for the dissimilarities between the clones. 


We realize that only a fraction of the total D. octaedra genome has been explored but in any case 
electrophoresis is an effective method to estimate genetic differences in animals that cannot be 
interbred. However, the kinds of enzyme systems explored may have an influence upon the 
interpreting data received. In the present study for instance, if one screens several monomorphic 
enzyme systems together with only Mdh and Lap ones, one gets 8clones with at least 
15 individuals, a few clones with some and three clones with only one worm in each. Introducing 
Idh and Pgi increases the number of unique clones to 18 and causes some scattering among the few 
worm clones but still leaves 5 clones with at least 15 individuals. Introducing Pgm and Est-2 
increases the number of unique clones to 80, the class of two-worm clones becomes as high as 30 but 
there still are five clones with at least 15 worms. The apparent conclusion is that introducing 
sensitive new systems will increase the unique clone category but still leaves the major groups 
intact. It is worth recalling that in table 4 we demonstrated that the different enzyme systems 
markers have significantly and also dissimilarly contributed to the delimitation of the clones. 
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Fig. 12. Efficiency of four polymorphic enzyme systems in delimiting groups of clones that respond to a given 
gradient a. In enzyme system A variants AI—A3 do not indicate any response for the clones with respect to the 
gradient. In enzyme systems B and C variants B1 and B2 as well as C1 and C2 indicate clear-cut ranges for the clone 
cathegories with respect to the gradient a. In enzyme D the three variants show rather trends than clear-cut differences 
between the clone groups. For further information see the text. 


There is, of course, an element of subjectivity in assessing electrophoretic material derived from 
animals that cannot be crossed and which, in addition are polyploid. In any case, the interpretation 
used here avoids lumping of otherwise similar-looking enzyme phenotypes in every case where a 
repeatable difference was seen. 


4.3. Clonal diversity in D. octaedra and other parthenogenetic earthworms 


The other electrophoretic surveys made on parthenogenetic earthworms concern Octolasion 
tyrtaeum [JAENIKE et al., 1980; JAENIKE & SELANDER, 1979, 1985], Dendrodrilus rubidus 
[JAENIKE et al., 1982], Eisenia nordenskioeldi [MALININA & PEREL', 1984] and triploid strains of 
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Fig. 13. Numbers of clones in the parthenogenetic earthworms Dendrobaena octaedra (based on the material in the 
present study) Dendrodrilus rubidus (data by JAENIKE er al. 1982) and Octolasion tyrtaeum JAENIKE ET AL. 1980) in 
rarefied samples according to the rarefaction method (see the text). Vertical bars correspond to the standard deviations 
(xi: S.D.). 


Allolobophora caliginosa [COBOLLI SBORDONI] et al., 1987]. The first two mentioned species can be 
directly compared with D. octaedra. 

O. tyrtaeum reproduces through obligate and D. rubidus through facultative parthenogenesis 
(MULDAL, 1952; OMODEO, 1952a, 1955). Clone diversity is much larger in D. octaedra than in the 
two other species (fig. 13) though based on fewer polymorphic enzyme systems explored. 

The relative proportion of clones with ten or more worms is 12.9% of all the clones and 29.5% 
of individuals in D. rubidus, 25.0% of all clones and 87.5% of individuals in O. ryrtaeum and in 
D. octaedra the corresponding figures are 5.4% for clones and 37.9% for individuals. 

As to unique clones, they comprise 0.1 96 of individuals in O. tyrtaeum and 12.5% of the total 
number of clones (JAENIKE et al., 1980). In D. rubidus the corresponding figures are 11.5% and 
51.596. The latter numbers are very close to the figures obtained here for D. octaedra, namely 
18.7% of individuals and 54.4% of clones. The data on D. rubidus are based on 5, those of 
D. octaedra on 6 and those of O. ryrtaeum on 10 polymorphic enzyme systems. Comparisons of 
clone diversities between triploid A. caliginosa [COBOLLI SBORDONI ef al., 1987] and D. rubidus 
[JAENIKE et al., 1982] populations account for a higher degree of genetic diversity among the 
former species. 

Both D. rubidus and O. tyrtaeum are introduced in North America, where they were studied. 
Thus the populations are relatively young and low clonal diversity may be due to a simple founder 
effect. Some of the types may have originated in America after the introduction; the facultative 
parthenogenesis of D. rubidus allows, at least in theory, recombination. 
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4.4. Reproduction and evolutionary aspects in D. octaedra 


JANDA & GavriLov (1939) concluded that D. octaedra does not reproduce through self- 
fertilization. On the basis of histological (CERNOsvITOV, 1930) and cytological evidence (OMODEO, 
1955) D. octaedra is regarded as a strictly parthenogenetic species. 

We have shown an abundance of D. octaedra clones in eastern Fennoscandia. There were only 
3 localities where all the worms within a sample were identical plus 7 others where only two overall 
genotypes were found within a sample of 5 worms. One might conclude that at least these localities 
are inhabited by clones. Any recombination event within a lineage would lead to decreasing levels 
of heterozygosity, which was not observed. No diploid populations have ever been recorded 
(OMODEO, personal communication) but based on morphological characters GATES (1974) 
hypothesized that they may occur in the Russian taiga belt. 

The D. octaedra worms are probably hexaploid or aneuploid with chromosome numbers close 
to hexaploidy (HONGELL & TERHIVUO, 1989) a level of ploidy common in very many populations 
(OMoDEO, 1955, 1957; CASELLATO & RODIGHIERO, 1972). Some worms may have all stages of 
spermatogenesis but still reproduce through apomictic parthenogenesis while others may show 
reduction of secondary male reproductive organs thus preventing amphigonic reproduction 
(TERHIVUO, 1988a; HONGELL & TERHIVUO, 1989). It is noteworthy that in many obligatorily 
parthenogenetic species there are rare males who arise from a variety of causes but they seem to have 
no reproductive function (e.g. SAURA, 1983). 

OMoDEO (1953) used the geographical distribution of D. octaedra (Greenland, different arctic 
islands, Europe and along the Mediterranean) as an evidence for the great antiquity of polyploidy 
and parthenogenesis in this morphologically variable taxon. He argued that we are confronted with a 
cluster of forms comparable to some complexes of apomictic plants. He used plate tectonics to 
determine the age of the disjunct populations and placed their origin in Tertiary. 

D. octaedra may have spread through passive transport over much of its present distribution 
area and the arguments of OMODEO can be debated; anyway the taxon is a versatile and eurytopic 
cluster of parthenogenetic clones. In the present study clone A is either an old form, a successful 
overall genotype or both. We do not need to affix any weight to its genetic constitution; the enzymes 
studied here serve only as markers. We have shown that it is a good "general purpose genotype" and 
may have given rise to a cluster of other forms through mutations. 

JAENIKE & SELANDER (1979) suggested that the two abundant clones of O. tyrtaeum represent 
"general purpose" genotypes but failed to show any ecological dissimilarities between them. They 
outlined the origin of parthenogenesis and polyploidy in a scenario where a dominant mutation first 
gives rise to parthenogenesis. Polyploidy results from chance hybridizations of unreduced 
parthenogenetic eggs. 

Polyploidy may also have come first: parthenogenesis offers a polyploid an escape from reduced 
fertility in polyploid meiosis. D. octaedra has apomictic parthenogenesis with only mitotic germ 
cell divisions. This is somewhat exceptional, as other parthenogenetic lumbricids have automixis 
(e.g. MULDAL, 1952; OMODEO 1952a and SUOMALAINEN et al., 1987). 

JAENIKE & SELANDER (1979) consider O. tyrtaeum to be an r-selected species that occupies 
patchy temporary habitats. In fact, these clones are succesful European colonizers of North 
America. As for D. octaedra it has the broadest habitat spectrum of all North European lumbricids 
(TERHIVUO & VaLoviRTA, 1978) and it has also highly successfully colonized many parts of the 
world (SrÓóP-Bowtrz, 1969). Understanding the genetic structure of a parthenogenetic taxon helps 
to overcome the old stereotyped notions of these forms as evolutionarily inferior "weeds". 
D. octaedra may possess signs of reduction in some characteristics but their patterns of clonal 
variation are organized in an interesting way that also seems to make good ecological and 
evolutionary sense. 


4.5. Coexistence of clones 


Electrophoretic studies on different parthenogenetic earthworm species have uniformly shown 
that several clones can coexist in the same place. JAENIKE et al. (1980) and JAENIKE & SELANDER 
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(1985) did not find any ecological or geographic differentiation between the two most abundant 
clones that, nevertheless, showed extensive temporal variability at a site. They considered three 
alternatives. One clone may be systematically replacing the other one in North America but the 
competitive exclusion has not run its course. They did not, however, find any clinal pattern in 
support of this contention. There was some slight indication of ecological differentiation, hence the 
localities in which clone A predominated tended to be moister than the others. The authors favoured 
clonal drift as an explanation. The two clones showed, however, regular co-existence over a wide 
area. 

We have demonstrated geographic clines and non-random distributions of clones over different 
biotopes and vegetation zones. The results agree with those of CHRISTENSEN (1980) and 
CHRISTENSEN ef al. (1976, 1977, 1978) who have shown clear-cut patterns in parthenogenetic 
enchytraeids; the results obtained here also approach the patterns observed in the polyploid and 
parthenogenetic populations of the weevil Otiorhynchus scaber [SAURA et al., 1976]. 

We have not attempted to study the composition of clones within a small area. A sample of a few 
hundred worms from a single locality would tell more about any partitioning of the habitat between 
different clones as well as about any temporal changes, given that it would be repeated at regular 
intervals. A worthwhile project would also be to study the composition of Octolasion tyrtaeum and 
Dendrodrilus rubidus populations in Europe. We predict that there would be far more clones present 
than in North America. 


4.6. Dispersal of D. octaedra to eastern Fennoscandia 


D. octaedra can be found in all parts of northern Europe (SróP-Bowirrz, 1969: TERHIVUO, 
1988b). Following the retreat of the Scandinavian ice sheet D. octaedra probably first colonized 
eastern parts of Finland from the U.S.S.R. The distribution of the clone group 4 indicates ties to the 
west. Dispersal to the Aland islands from Sweden and to central Finland from northern parts of 
Sweden could explain this tie. 

D. octaedra can occupy even small skerries in the Áland archipelago. Cocoons or individuals 
stuck to leaves or pieces of wood debris may have been carried by waves all over the sea. This holds 
true at least in the brackish water of the Baltic (HAEGGSTRÓM & TERHIVUO, 1979). Sea currents 
running along eastern parts of the Baltic may have facilitated its dispersal over the Gulf of Finland 
from Estonia and nearby areas. Evidence for this kind of dispersal is available for many insect 
species of the area (PALMÉN, 1944). 

Agricultural tools and machines can effectively spread earthworms and their cocoons. This 
explains some but not all of the dispersal over areas with agricultural practices. Parthenogenetic 
forms are efficient dispersers, as a viable cocoon suffices to establish a new population. D. octaedra 
individuals, probably more often cocoons, may face many kinds of environmental factors 
combinations (frost, low pH, raw humus soils, wet or dry soils etc., see TERHIVUO, 1989) 
adversities related to long distance dispersal jumps. Much of the variation observed in natural 
populations can be explained through migration patterns while the remainder probably can be 
accounted for by local differentiation of populations through mutations. 


5. Summary 


D. octaedra is a polyploid earthworm species that reproduces through apomictic parthenogenesis. Samples 
consisting of about 5 worms each were collected from 85 localities situated 20—40 km apart along two transects, one 
of which extends for over 1,000 km from north to south and the other over 600 km from east to west across southern 
Finland. A total of 428 worms were subjected to enzyme electrophoresis. Six loci proved to be variable in the material. 
D. octaedra showed wider clonal variability than observed in other parthenogenetic earthworms studied up till now. A 
total of 147 different clones were found and 80 of these were unique, represented by one worm only. The common 
clones and clone clusters (defined by the number their enzyme types deviating from those in the most abundant clone) 
showed distribution patterns that can be correlated with the biotopes and vegetation zones along the N-S transect. 
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Clonal diversity decreased slightly towards the north in rich deciduous forests and in dry coniferous heath forests but in 
meadows it increased. Areal as well as between-biotope similarities of clone pools were low indicating the high 
evolutionary potential of the species even in marginal parts of its geographical range. Evidence for post-glacial 
migration patterns influencing genetic diversity of D. octaedra populations in different parts of eastern Fennoscandia 
is discussed. 
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Synopsis: Original scientific paper 
TERHivUO, J., & A. Saura, 1990. Allozyme variation in parthenogenetic Dendrobaena octaedra (Saw.) 
(Oligochaeta: Lumbricidae) populations of eastern Fennoscandia. Pedobiologia 34. 113—139. 

Genetic variability in marginal parthenogenetic Dendrobaena octaedra populations in eastern Fennoscandia was 
explored electrophoretically. Six polymorphic loci were used to determine the overall genotype for the 428 
D. octaedra individuals sampled in different biotopes and parts of the study area. D. octaedra showed wider allozymic 
variability than the other parthenogenetic earthworms studied up till now. The frequent clones and clone groups 
showed unequal distribution patterns in relation to biotopes and regions. Clone diversity decreased only slightly 
towards the north and clone pools were very dissimilar both between the biotopes and between the areas indicating the 
great evolutionary potential of the species in marginal areas where dispersed postglacially. Clues to the migration 
pattern of the clones to the area were also found. 

Key words: Earthworm, parthenogenesis, electrophoresis, genetic variation, enzyme. clone. 
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